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ABSTRACT 

I discuss the status of ongoing work to determine numerically the elastic ttt: scattering 
amplitude at order p 6 in the framework of chiral perturbation theory. 

1. Introduction 

The interplay between theoretical and experimental aspects of elastic irn scattering 
is illustrated in figure 1. On the theoretical side, Weinberg's calculation! of the scattering 
amplitude at leading order in the low-energy expansion gives for the isospin zero S-wave 
scattering length the value o/Tq = 0.16 in units of the charged pion mass. This differs from 
the experimentally determined valueH a{] = 0.26 ± 0.05 by two standard deviations. The 
one-loop calculation! enhances the leading order term to a[j = 0.20 ±0.01 - the correction 
goes in the right direction, but the result is still on the low side as far as the present exper- 
imental value is concerned. To decide about agreement/disagreement between theory and 
experiment, one should i) evaluate the scattering lengths in the theoretical framework at 
order p 6 , and ii) determine them more precisely experimentally. My talk was concerned 
with the the former issue, indicated by the double arrow in the left column of the figure. 

On the experimental side, several attempts are underway to improve our knowledge 
of the threshold parameters. The most promising ones among them are semileptonic 
decays with improved statistics (E 865iand KLOE&), and the measurement of the pionium 
lifetime (DIRAOI) that will allow one to directly determine the combination | a§ — a§ | of 
S'-wave scattering lengths. It was one of the aims of the workshop to discuss the precise 
relation between the lifetime of the pionium atom and the tttt scattering lengths - I 
refer the interested reader to the numerous contributions to this workshop for details. In 
addition, Pocanic! has presented a nice review of the determination of scattering lengths 
from hadronic processes. 

*' Talk given at the International Workshop "Hadronic Atoms and Positronium in the Standard Model" , Dubna, 
Russia, May 26-31, 1998. To appear in the Proceedings. 
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As for theory, I will describe in the following sections our current effort^ to deter- 
mine the scattering amplitude numerically at order p 6 in the framework of chiral pertur- 
bation theory (CHPT). As the title of my talk indicates, this work is still in progress - 
therefore. I wil eive an overview of what is going, omitting any details. For related work, 
seeRefsSi'M. 
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Figure 1: Progress in the determination of the elastic tttt scattering amplitude. References 
are provided in the text. 



2. Chiral perturbation theory 

To set notation, I briefly describe the framework use to evaluate the tttc scat- 
tering amplitude. I consider the isospin symmetry limit where the two lightest quark 
masses are equal, m u = m<2 ^ 0, and where the electromagnetic coupling is set to zero, 
O-qed = 0. I refer to Miiller'slil contribution for a discussion of isospin violating effects 
in hadronic amplitudes. The effective lagrangian that describes the interaction of pions 
is given by a string of terms, £ e g- = C 2 + ^£4 + ti 2 ^ + • - • , where C n contains mi 
derivatives of the pion fields and m <> q uark mass matrices, with mi ± 2rri2 = n (here, I 
consider the standard counting rulesli-i' see section 7 for a generalization thereolc3). 
The low-energy expansion corresponds to an expansion of the scattering matrix elements 
in powers of H. 

In the following, the low-energy constants (LEC's) hidden in the effective la- 
grangians C n play an important role. In £ 2) there are two of them, the pion decay con- 
stant in the chiral limit (F) and the parameter B, which are related to the condensate by 
F 2 B = — (0\uu\0). In the loop expansion, these two parameters can be expressed in terms 
of the physical pion decay constant F n ~ 92.4 MeV and of the pion mass, M T = 139.57 
MeV. The scattering amplitude contains, in the two-loop approximation, in addition 
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several LEC's ocurring in £ 4 and in £ 6 , 



£2 



h, h, hi h 
fx,..., fa 



> occur in tttt — > tttt (two-loop approximation). (1) 



These LEC's are not determined by chiral symmetry - they are, however, in principle 
calculable in QCdH. 

3. tttt -> tttt in CHPT 

Elastic tttt scattering is described by a single Lorentz invariant amplitude A(s, i, u), 
that depends on the standard Mandelstam variables s,t, u (I use the same notation as 
SainiocM- I refer the reader to his contribution for details). Here, we are concerned with 
the loop expansion of the scattering amplitude, 



A( S ,t,u)={ A > + ^ + * + °< p8 » ' (2) 
tree 1 loop 2 loops 

where A n is of order p n . In the following, I denote with the symbol A x the first three 
terms in (H), 

A x = A 2 + A 4 + A 6 . (3) 

A dispersive evaluation of A x has been performed in Ref.HH. That calculation is not 
sufficient for the present purpose, because several tadpole diagrams that occur in the 
loop expansion have not been evaluated in that work. What is needed for the analysis 
outlined below is the complete two-loop expression of A x . The relevant calculation has 
been performed in Ref.13. The explicit form of the amplitude is not needed here, and I 
refer the interested reader to Sainio's contribution!^! for explicit formulae. 

For the following discussion, it is useful to bring the amplitude into a canonical 
form. First, one makes use of the fact that A x can be expressed^ in terms of three 
functions Wo,i,2> that depend on a single kinematical variable (s,t or u) and have only a 
right hand cut. Their imaginary part is given by the imaginary parts of the three lowest 
partial waves if, e.g., ImW/o = Imi[j, etc. We now write a four times subtracted dispersion 
relation for W , 

W = a + ai s + a 2 s 2 + a 3 s 3 + — / — . (4) 

7T J4,M% X 4 X — S 

and similarly for Wi^. Inserting these representations into the expression for A x gives 

AX — pX 1 AX 

P x = (3 X + {3 x s + f3*s 2 + f3l{t-u) 2 + f3*s z + f3*s{t-u) 2 . (5) 
The coefficients 0f contain the low-energy constants (0), 

Pn = P%(h> ■ ■ ■ Jtin, ■ ■ ■ ,r 6 ). (6) 
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Here and in the following, I drop all dependence on F n and whose numerical value 
is taken from experiment. It is the aim of the enterprise to work out the consequences 
of the requirement that A x is a good approximation of the true amplitude at low energy. 
This requirement allows one in particular to work out the two S-wave scattering lengths 
a,Q and o,q, provided that we can determine the LEC's ocurring in the chiral amplitude. 



4. Effective couplings 

The chiral amplitude A x contains two sets of LEC's: 



set 1 In this set we put the four LEC's 



set 2 



hi hi ?"5j r 6 



(7) 



that are related to the momentum and scattering angle dependence of the amplitude: 
the polynomial P x contains terms of the form l\s(t — u) 2 , l%{t — u) 2 , fss 3 , fes(t — u) 2 . 
As a result of this, these couplings can be determined from data above ~ 800 MeV, 
see below. 

This set assembles the LEC's whose effect disappears in the chiral limit, because 
they are multiplied with the square of the pion mass. Examples are the terms 
P x . As nature does not allow us to vary the pion mass in the 
laboratory, we need input from outside tttt scattering to determine these couplings. 
I denote these with the symbol B, 
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^3,4j r l,2,3,4 



Using large- N c arguments, the constant I3 can be obtained from the chiral expansion 
of the ratio M^/M 2 0, whereas I4 is related to the scalar radius of the pionS H. 
Finally, the effect of ^2,3,4 is suppressed by powers of the pion mass - a rough 
estimate like the one provided in Ref£H (based on resonance saturation) thus suffices. 



5. Matching 

From the above discussion, we find that the coefficients f3 x in the polynomial P x 
of the chiral amplitude (0) are of the functional form 

ffi = l3X(h,l 2 ,f 5 ,f 6 -,e). (9) 

In the following, I consider the couplings G as given. 

Turning now to experiment, I write the corresponding amplitude A exp as 

A ex V = pexp + A exv ( 1Q ) 

where P exp has the same functional form as P x , with (3 X —>■ (5 exp . As we will see in the 
following section, the Roy equations^! fix these coefficients in terms of data above ~ 800 
MeV and of the scattering lengths a[], a^, 

P exp ~P e n xp (a° Q ,a 2 ). (11) 



4 



Furthermore, the amplitude A^f coincides with Af nt up to terms of order p 8 . We now 
require that the chiral amplitude agrees with the experimental one near threshold , as a 
result of which one obtains the relation 



P x (s,t,u) = P exp (s,t,u) 
matching condition 



'12) 



In terms of the coefficients f3 n , this matching condition amounts to the six relations 

ffiih, l 2 , f 5) f 6 ; 6) = (3 e n xp (a° 0} ojj) ; n = 1, . . . , 6 (13) 
for six unknowns. Solving for these, we can determine the quantities 



a8,ag;Zi,Z 2 ,r 5 ,r B . (14) 



The remaining threshold parameters may be obtained from the Wanders sum rulesl!> M. 

6. Roy equations 

It remains to show that indeed the polynomial coefficients /3^ xp are determined in 
terms of the two S-wave scattering lengths a[j, and of data above 800 MeV. 

For this purpose, one consider the partial waves tj. As shown by RoyH, their real 
and imaginary parts are related through 

2 00 poo 

Retj = c I l (s) + J2Ef , teKw («, *) Im 4 (*) (15) 



Ii=0l>=0 MM % 

in the region — 28M2 < s < 60M2. Here, the quantities K(f are known kernels, whereas 
the subtraction polynomial c\ contains the two scattering lengths the only free 

parameters. Note that the relation (|T5|) is linear in the partial wave amplitudes - it 
is therefore also true in ordinary perturbation theory, order by order in the coupling 
constant (modulo subtractions). Next, we observe that, at low energy, the most important 
contribution to the imaginary parts stems from the S- and P-waves, because the higher 
waves are suppressed. Imtf = 0(p 8 ) ; I — 2, 3, . . . . For this reason, it is useful to perform 
the splittingpJ 

A exp = A SP + A R , (16) 

where the crossing symmetric amplitude A$p contains the exact S— and P-wave absorp- 
tive parts, whereas Ap has only absorptive parts from I = 2, 3, . . . waves. Therefore, at 
low energies, the latter contribution can be approximated by a polynomial in s, t and u. 
This allows one finally to write A exp in the form (|10D , where P exp is determined from Ap 
and through the imaginary parts of the three lowest partial waves. I now show how one 
may determine these imaginary parts from experimental information above 800 MeV. 

Using ([IB]) and expanding A R as described, the partial wave relations ( |T5|) become 
for the S- and P waves 

Ret m = c m + V F dxK mn (s, x)lmt n (x) + d m (s) ; m = 0, 1, 2 (17) 
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where 

(t°, t{, tl) -> (t , *i, t 2 ) ; (eg, cj, eg) -> (c , ci, c 2 ) . (18) 
The kernels are Kqo — ^oo , etc. Finally, the so-called driving terms d m are obtained from 



the amplitude Ar. The crucial stepc3 is to invoke unitarity at this stage, 

tm = ^ (i-4M2 /s) i/ 2 (Vmexp2t5 m - 1) . (19) 

We now assume that 

• the elasticities r] m are known for all s 

• the 5 m are given above a cutoff energy E cut ~ 800 MeV 

• the driving terms are known below E cut . 

At given a[j, Oq, the equations fll7D then become a system of nonlinear integral equations 
for the unknown 5 m in the interval 2M n < y/s < E cut . This system is known as 'Roy 
equations' in the literature. Solving theseQ, one obtains the S, P partial waves at low 
energies and thus the polynomial P ex P^ as promised. 

7. Generalized chiral perturbation theory 

The above framework assumes the standard scenario of chiral counting, where 
the condensate is the leading order parameter, considered a quantity of order one. In 
recent years, an alternative picturelfhas been implemented^! in an effective lagrangian 
framework. In that picture, the condensate may be small or even vanishing. Many 
predictions are lost in this prominent one being the Gell-Mann-Okubo mass 

relation, that does not have a natural explanation in this framework. Furthermore, the 
scattering length a° cannot be predicted either - it is a free parameter, related to the size 
of the condensate. The philosophy of this approach is to let experiment tell the size of 
a[}, Oq, which allows one then to determine the size of the condensate. I refer the interested 



reader to Stern's talk in MainzEa for further details and for references. 

It is not easy to distinguish phenomenologically a small condensate from the stan- 
dard case with the present precision of low-energy experiments". It is, however, expected 
that new precise measurements of the irn amplitude at low energy will shed more light 
on the issue. 

On the theoretical side, a recent interesting investigation by Knecht and de RafaelH 
has shown that, in vector-like gauge theories like QCD, and at large values of N c , the 
ordering pattern of narrow vector and axial-vector states is correlated with the size of 
possible local order parameters of chiral symmetry breaking. The authors find that, 
from a duality point of view, the option of a vanishing condensate seems unnatural. I 
refer the reader to their work for more details. Recent lattice calculations^' do not 
support a small condensate either. Indeed, Ecker^l has compared the results of Ref.^3 
on the quark mass dependence of the meson masses with the predictions of standard 
and generalized CHPT and concludes that ". . . lattice QCD is incompatible with a small 
quark condensate". The authors of Ref.iH have made a considerable effort to pin down 
systematic uncertainties in a direct evaluation of (0|gg|0) on the lattice in the quenched 
approximation. They come up with the standard value, with remarkably small error bars. 

^For E cut ~ 800 MeV, the solution is uniquely. 



6 



8. Summary 



Status 



The analytic form of the two-loop amplitude in standard CHPT is knownS. 

We can numerically construct the solution of the system (f[7l) for given scattering 

lengths, input data and driving termsli3. 



The determination of the uncertainties in the final values of the parameters ([14]), in 
particular in the scatterin g le n g ths, is in progress^. 
For related work, see Refs 0' S © — . 



Outlook 



gchptH provides a framework where the condensate may be small or even van- 
ishing. 

In this scenario, the isospin zero S'-wave scattering length a[j can be large, in contrast 
to the standard case presented here, and in contrast to lattice calculations© Hi 11 
in the quenched approximation. 

We all hope that DIRAci, E865I, KLOfJ and CHAOsH will soon provide addi- 
tional information on the issue from the experimental side. 
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